Motion of Massive and Massless Test particles 
in Dyadosphere Geometry 

B. Raychaudhuri* F. Rahamanj M. Kalam-'-and A. Ghosh^ 



Abstract 

Motion of massive and massless test particle in equilibrium and non-equilibrium 
case is discussed in a dyadosphere geometry through Hamilton- Jacobi method. 
Geodesies of particles are discussed through Lagrangian method too. Scalar wave 
equation for massless particle is analyzed to show the absence of superradiance. 

1 Introduction 

It is now well established that a black hole is uniquely characterized by its mass M, 
charge Q and angular momentum J. The exterior solution for a rotating charged black 
hole is given the Kerr-Newman metric while that for the uncharged one is given by 
Kerr geometry. Further, if the black hole is charged but non-rotating, then the solution 
is Reissner-Nordstrom metric, while for uncharged case the black hole is represented 
by the Schwarzschild geometry. 

Recently Ruffini [Ij and Preparata, Ruffini and Xue [2j discussed a special region 
just outside the horizon of charged black holes where the electric field goes beyond 
its classical limit. This is a situation where effects of vacuum fluctuations should be 
considered. When the electromagnetic field exceeds the Heisenber-Euler [3J critical 
value for e^e~ pair production, in a very short time of the order ~ 0{h/mcP), a very 
large number of e^e~ pair is created there. This reaches in thermodynamic equilibrium 
with a photon gas. The period of overcritical state is of approximately 10'' s. [2]. The 
dissipation occurs through a process involving the dyadosphere of a black hole. 

Taking into account the one-loop correction due to QED to the first order of ap- 
proximation the space time geometry has the form 

dr^ 

ds^ = -/(r)df2 + + {de + sin^ Odcj)'^) , (1) 

where 

/(O = 1 - — + 7^ - (2) 

o" is a coupling constant. The metric reduces to usual charged non-rotating RN metric 
in the limit a = 0. 
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Preparata, Ruffini and Xue [2] discussed some inequalities which must be satisfied 
during the forming process of a dyadosphere of a black hole. The authors put forward 
a theoretical model on the concept of dyadosphere explaining Gamma- Ray-Bursts phe- 
nomenon. Lorenci et. al. [Ij discussed the expression for the bending of light in the 
dyadosphere of the charged black hole. 

In this paper we study the motion of test particle in the field of the dyadosphere 
of a black hole. The procedure outlined by Rahaman [5] is followed here. In Sec. [2] 
we discuss both static and nonstatic cases for charged and uncharged test particle [3]. 
Sec. [3] will discuss the motion of test particles in static equilibrium while Sec. [H will 
describe the motion of test particles in non equilibrium state. We discuss the wave 
equation for a massless scalar wave in Sec. [5] to split and solve the radial and angular 
equations. Absence of superradiance is expectedly observed in this case. Finally in 
Sec. El geodesies of dyadosphere metric using Lagrangian technique is discussed. 



2 Trajectory of Test Particle through Hamilton- Jacobi 
Formalism 

We first consider a test particle with mass m and charge e moving in the modified 
gravitational field of a RN black hole due to the presence of a dyadosphere (henceforth 
called dyadosphere field). The Hamilton Jacobi equation for the test particle is given 
by [6] 

where gik is the background geometry, i.e. the metric ([I]) and 

Ai = > gauge potential. (4j 

r 

S is the usual Hamilton's characteristic function. 

For the metric ([T]) the explicit form of the Hamilton-Jacobi equation is given by [5] 

To solve this partial differential equation, an ansatz is proposed [5]. 

S{t, r, 9, cj)) = -E.t + 5i(r) + 52(0) + J.cP. (6) 

This will enable us to separate the variables. E and J are identified as the energy 
and angular momentum of the particle. The substitution of this ansatz in the previous 
equation yields 

dS___^ dS _ dSi 

dt ' dr dr ' ,„■, 

dS_ _ dS2 dS_ _ ^'^ 
'de ~ ~de' dQ ~ ■ 
Substituting we obtain the following equations 
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With all these, we find the following solutions for the functions Si{r) and S2{9) as 



Si{r)=e 



jp eQ 



P 



2 2 

m p 

T7^ 



1/2 



dr, 



(12) 



and 



S2{e) = e I [p^-J^csc^efdO, 



(13) 



where e = ±1. 

To determine the trajectory of the particle through Hamiltonian Jacobi method, we 
consider [5] 

dS dS dS 

— — = constant, — — = constant, — — = constant. (14) 
oE oJ oP 

Without any loss of generality the constants may be considered to be zero. Then, 
from Eq. ^ we obtain 

-t (15) 



dE 



which gives 



t = el'(E-^^ 



f fr^ 

€ j (J csc^ 6) [p^ - J2 csc^ e] ^'"^ de, 



-1/2 



dr. 
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The Radial velocity is particle is obtained from Eq. (fT6 



(16) 
(17) 

(18) 



d'' _ f2 fjp eQ\ ^ 



— (e- ^ -— 



/ 



(19) 
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dr 

The turning point of the particle is — = and the potential curves are 

at 



E- 



eQ 



(20) 



Solving 



m + "2 



The effective potential is given by 



m mr • - - ^ m?r'^ 



1/2 



as an explicitly function of r only. A representative curve is given below (Fig. [T]). We 
have considered the specific charge of the particle e/m = 1, and p/m = 1.5 and a = 0.5 
for convenience. This is a plot for a extremal case where Q = M. We set them to unity. 




Figure 1: V{r) for e/m = 1, Q = AI = 1, a = 0.5 and p/m = 1.5 

Different curves may be drawn depending upon different values of the parameter. 
We give another plot where M = 1.5 but Q = 0.5 (Fig. [2]). 

In a stationary system, E and so V{r) must have an extremal value. Hence the 
value of r for which the energy attains its extremal value is given by the solution of 

dr 



This gives 



mr^ 2y/f \ m'^r'^ J dr \ m^r"^ J m^r^ 
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Figure 2: V{r) for e/m = 1, Q = 0.5, M = 1.5, a = 0.5 and p/m = 1.5 



We obtain 



1/2 



2 V m^r^ / (?r 



5/ ^ F 



With the expression of /(r) the equation becomes 



eQ ( 2M_ Q2 

1 ~t~ o 



m 



4\ 1/2 



5^6 



1 + 



1/2 



r 5 / V m^r^ 



This expression is used in further discussions below. 



3 Test particle in Static Equilibrium 



2M g2 



5^6 



(22) 



In static equilibrium, momentum p in Eq. (|22|) must be zero. So, the value of r for 
which potential wih be an extremal is given by 



m 



2M Q2 ^q4x1/2 



5^6 



r 5 r"' 



Simphfying and arranging 
^^Q^\ „io 



2MQ' I 1 - 4 V' + fl - 4 1 Q'^' 
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This is a 10th order equation in r. If — ^ > — then the last term is negative, 

m Q 

(? 

i.e at least one positive root exists. If on the other hand, — ^ — — 2"' then the above 

m Q 

equation reduces to a ninth degree equation with negative last term implying that a real 
positive root exists. In other words, the particle will be trapped by the dyadosphere 
field. 

Note that if the particle is of unit specific charge, e/m = 1 then the equation 
becomes 



(24) 



i.e. for an extremal RN field, the equation becomes a 5th degree equation on r. The 
function is plotted below (Fig. [3]). Again we demonstrate the nature of the function 
with M = 1, e/m = 0.5, Q = 0.7 and a = 0.5. The function is monotonically increasing 
outside the range specified. 
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Figure 3: M = 1, e/m = 0.5, Q = 0.7 and a = 0.5 



4 Test particle in Non equilibrium state, p 
4.1 Case I — Uncharged test particle (e = 0) 

We first discuss the case of uncharged test particle e = 0. In the nonequilibrium case, 
p ^ 0. The expression (|22p becomes 



3aQ4W \ p' ( 2M Q2 ^g4x 
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Simplifying and rearranging we obtain an algebraic equation of degree 7. 

Mr^ - f q2 ^ 4 V' + 34^^' - '^^Q'^' + I'^Q'^' + 7^Q'4 = (26) 
\ m^y mb2 5 5 

Since there are four changes of sign in the above equation, then by Descartes rule 
of signs, the equation can have 4 positive roots. 

4.2 Case II, Charged particle (e 7^ 0) 

In case of charged particle e 7^ and the Eq. (j22p becomes 



r 



(271 

_ f2 A _ ^ + «! _ ^«!V'Vi + ^V'^ 0. 

This is an equation of degree 14 in r. When explicitly written this expression 
contains 60 terms. In terms of descending power of r it reduces to 



M- - _ 2M f g2 + 4 - ^ 1 



+ 4 f 2M2 + 4M + - 2Q2 + 4 - (28) 



... + Oxr + |^cj2g8 = o 
5 

In this case too, if the specific charge of the test particle e/m is greater than M/Q 

then the last term is negative, i.e. at least one positive root exists. If on the other hand 

e M 

— = -771 then the above equation reduces to a 13 degree equation, with negative last 
m Q 

term if — 9 ^ ^ implying the existence of a real positive root. In these cases 

m 

the particle will be trapped by the dyadosphere. 

5 Solution of massless Scalar Wave Equation in Dyado- 
sphere metric 

In this section we shall analysis the scalar wave equation for dyadosphere geometry fol- 
lowing Brill et. al f9] who first showed that variables can be separated for Kerr geometry. 
(We consider massless scalar field here; extension to massive field is straightforward and 
devoid of any further conclusion.) The wave equation for a massless particle is given 

.-"^^(."V"^)x^O (28) 

Here gf^u is given by Eq. ([1]). If we now put all the values we obtain the following 
equation 
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Where A = - 2Mr + - 

This equation be solved by separation of variable with the ansatz 

X = e-''^'e'""''R{r)e{e) (31) 
Putting this in the wave equation we obtain 

r^sin6' 2 si^^ ^ /. dR\ 19/. „dQ\ m? 



Thus the variables can be easily separated. The radial equation reduces to 

a|: f a|^1 + (rV - AA)i? = (32) 



dr V dr 
The angular part becomes 

1 d f . BQ\ 



m 



2 



. . sine— --^e + Ae = (33) 

sm 6 o9 \ od J sur 9 

With a substitution x = cos 9, the equation becomes 

(i-^')3z^-2x3-: - ^-r^-A = (34) 



dx^ dx \1 

If we now write X = l{l + 1) and demand that / takes only integer values then the 
equation 

is Associated Legendre equation and the solution is give by the associated Legendre 
polynomial P[^{x) and is expressed as 

er(cose) = Prix) = (1 - ^T^' £^ - 1)' . (36) 

5.1 The Radial Equation: Absence of Superradiance 

Superradiance is the wave analogue for the Penrose process on Black Hole. If a bosonic 
or fermionic wave is incident upon a black hole, normally the reflected wave carries 
less energy than the incident wave. But under certain condition the transmitted wave, 
absorbed by the black hole carries negative energy into the black hole making the 
reflexion coefficient for the wave greater than unity. This means that the reflected wave 
will carry more energy than the incident wave. This phenomenon is called superradiance 
by Misner [8j and also analyzed by Zel'dovich and Starobinsky [TOlITT]. Through this 
process energy can be extracted from a black hole in expense of its angular momentum. 
The condition is given by [12] 



< u; < mnn, (37) 

where i^n is the angular velocity of the horizon ^ . Chandrasekhar |7j has analyzed the 
phenomenon for Kerr geometry in characteristic details and showed that this occurs 
only for incident waves of integral spins, i.e., for scalar, electromagnetic waves and 
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gravitational cases. He showed its absence for the fermionic waves, i.e. Dirac wave or 
neutrino waves. Basak and Majumdar discussed this phenomenon for acoustic analogue 
of Kerr Black Hole [13\ I14j. Clearly for a non-rotating black hole superradiance should 
be absent as we demonstrate precisely in this case of dyadosphere geometry. 
The radial equation is given by 

^Tr i^f) + - + '^^^ ^ = ^''^ 

Let us introduce the familiar r* coordinate (the tortoise coordinate) defined by 

dr* 

dT " A' 

thus giving 

A d 2 d 
A— = r^ — . 

dr dr* 

Note that though the variable r* is defined in the same manner as in Schwarzschild 
or in Kerr metric, in this case the variable is non-integrable. Still the basic purpose 
is satisfied, the coordinate spans over the real line and pushes the horizon to minus 
infinity. 

The introduction of another function u(r) = rR reduces the radial equation in a 
more familiar form 



4 - A 2A2 2AM , A ^ 

r^'^-—^ ^ + (2 + K^ + i));:4-^' 



Thus a potential barrier remains where 



0. (39) 



At horizon (A ^ 0), the radial equation becomes 

^ + .^u, = 0, (41) 

with V{r) = —ij?. 

Now asymptotically, r — > oo. The equation has the same form as in the previous 
case 

^ + -^-00 = 0. (42) 

Thus uu = Uoo, where uu is the radial solution at horizon and Uoo is the solution 
at cxD. This equality shows that for a dyadosphere metric there is no 'phenomenon of 
superradiance for an incident massless scalar field. This is expected as dyadosphere 
geometry is nonrotating. 



6 Analysis of Dyadosphere metric using Lagrangian Method 

The dyadosphere geometry evidently allows two timelike Killing vectors dt and d^j). 
Thus there must be two constants of motion. Note that the geodesic in space time can 
be derived from the usual Lagrangian \J\ 

dx* dx-^ , , , 

2£ = f«d7df 



9 



Here r is some affine parameter along the geodesic. For time-like geodesic r may 
be identified with the proper time s of the particle describing the geodesic. 
For the Dyadosphere metric the Lagrangian is thus given by 

2C = f(r)i^ - -^r - r^e^ - r^, sin^ 94>'^ (44) 

with the dot meaning the differention with respect to r. 
The canonical momenta are given by 

dC 

dqj 

The momenta in this case are given by [7] 

P. = f = fir)i, 

dr f{r) 



(45) 



= = r sm d(p. 

The Hamiltonian is given by 

n=pti- (prr + Pee + p^ij - C. (46) 
Putting the values of the conjugate momenta from Eq. (|45|) we find that 

n = c, 

which means the problem has no potential energy term. 

Rescaling of the affine parameter r allows us to set 2£ = 1 for time-like geodesies. 
For null geodesies £ = 0. 

As the metric allows two timelike Killing vectors dt and dff^ so we shall have two con- 
stants of integration, namely the conjugate momenta corresponding to the coordinates 
t and (f). Thus 

^ = — = n 
dr ~ dt ~ ' 

(47) 

dr dcj) 
Thus we find the first constant of motion 

Pt = /('")~i~ = constant = £'(say). (48) 
dr 

The metric is asymptotically fiat. So the quantity E can be identified to the energy at 
infinity. 
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The second constant of motion is given by 



P(f, = sin^ 6-^ = constant. (49) 
From the equation of motion we obtain 

Note that if we assign Q = -k jl when = then = 0. This means that the geodesic 
is described in an invariant plane which may be distinguished by = 7r/2. Thus we 
obtain 

odt/) 

= r — = constant = L (say). (51) 

Here L denotes the angular momentum about an axis normal to the invariant plane. 
Putting all these constants we obtain 

— = 2/: = +1 or 0. (52) 



fij) f{r) r 

depending whether it is a timelike or null geodesic. 
We find 



r^+{^ + 2C]f{r) = E\ (53) 



Thus we can identify the following quantity at the left hand side of the equation as 
effectively to a potential energy term fl\ pl02] 



r2 I ^ 



KE=[yI + ^^^y{r). (54) 
6.1 Radial timelike Geodesic, £ = +1, L = 

The effective potential for radial geodesic (zero angular momentum) is given by 

,55, 

which is equal to the original dyadosphere metric. The effective potential is plotted 
(Fig. g]) here for a black hole of M = 1.5, Q = 0.04, and a = 0.5. 



6.2 Null geodesic 

For null geodesic, £ = and the effective potential is given by 

^el = ^f{r). (56) 

Thus the effective potential for a photon with zero angular momentum vanishes. 
The effective potential is plotted (Fig. [5]) for a black hole of M = 1.5, Q = 0.7, cr = 5 
and L = 0.7. 
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Figure 4: V^s with M = 1.5, Q = 0.7 and a = 0.5 




Figure 5: V^s for L = 0.7, M = 1.5, Q = 0.7, and a = 0.5 

7 Discussion 

In this paper we discussed motion of test particle in the dyadosphere metric using 
Hamilton-Jacobi techniques. Some representative curves for effective potential is plot- 
ted. (All the curves are drawn by gnuplot plotting program). Several plots can be 

generated by changing the values of the parameters to see the nature of the function 
in different cases. We also discussed geodesies through Lagrangian method. The wave 
equation for a massless particle also is studied to show that the phenomenon of super- 
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radiance is absent in this case. 
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